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Abstract 

Events with four distinct jets from e+e~ collisions, collected by the OPAL detector at centre-of- 
mass energies between 130 and 184 GeV, are analysed for a peak in the sum of dijet masses. This 
search is motivated by the ALEPH Collaboration's observation of a clear excess of events with 
dijet mass sums close to 105 GeV in data taken at centre-of-mass energies of 130 and 136 GeV 
in 1995. We have observed no significant excess of four-jet events compared to the Standard 
Model expectation for any dijet mass sum at any energy. Our observation is inconsistent 
with the excess observed by ALEPH in 1995. Upper limits are determined on the production 
cross-section clS db function of the dijet mass sum. 
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1 Introduction 



In a run of LEP in 1995 at centre-of-mass energies of \fs = 130 and 136 GeV, the ALEPH 
Collaboration observed ||l]| an excess of events with four distinct jets compared with the Stan- 
dard Model expectation. Such an excess could be due to the production of new particles X 
and Y, each decaying into two hadronic jets in the process e+e~— >^XY— >^ four jets. The two 
particles could have equal or unequal masses. Grouping the jets into pairs, calculating their 
pair invariant masses Mij and Mm, and selecting the combination yielding the smallest mass 
difference AM = \Mij — Mki\, ALEPH observed a clustering of nine events in a mass window 
6.3 GeV wide centred around M = Mij + Mki ~ 105 GeV, with a Standard Model expectation 
of 0.8 events in this mass window. The choice of the combination with the minimum AM 
would tend to favour the selection of particles of equal mass or with a small mass difference. 

In response to this observation, the OPAL Collaboration performed an analysis that closely 
followed the selection of Reference [|l[. In its 130 and 136 GeV data from 1995, OPAL observed 
seven events with M between 60 and 130 GeV, with an expected Standard Model background of 
6.4±0.6 events. In the signal region indicated by ALEPH, OPAL observed one event, consistent 
with the expected Standard Model background of 0.8 ±0.2 events 0. The estimated efficiency 
of the OPAL analysis and the resolution on the dijet mass sum are similar to those obtained by 
ALEPH. Consequently, the OPAL detector would be expected to have a sensitivity comparable 
to the ALEPH detector for a four-jet signal should one exist. The ALEPH Collaboration also 
reported a slight excess at the higher centre-of-mass energies of 161 and 172 GeV 0. The 
DELPHI Collaboration observed no significant peak at 105 GeV in a similar analysis using 
their 1995 data at 130 and 136 GeV [Q. The L3 Collaboration likewise reported no excess of 
events in the indicated mass window for y/s = 130-172 GeV ||^. Nonetheless, there has been a 
great deal of theoretical speculation on the cause of the excess observed by ALEPH 0, |]. 

In 1997, LEP made short runs at y/s = 130 and 136 GeV with an integrated luminosity 
similar to that of 1995 at these centre-of-mass energies to test again the signal hypothesis. We 
add these data to our sample described in Reference [0], and also include data collected at 161, 
172 and 183 GeV. To search for the class of events observed by ALEPH in a model independent 
fashion, we have performed analyses on the OPAL data as close as possible to the ALEPH 
analyses at these energies P]. However, above the kinematic threshold for W-pair production, 
a veto is imposed to suppress this new source of background and results are presented with and 
without this requirement. 

The comparison of results of the OPAL emulation of the ALEPH selection to the ALEPH 
observation of an excess does not depend upon the underlying model of possible new physics 
if only the number of observed events is compared. In the context of a model of the process 
e"'"e~^XY^ four jets, a separate analysis is also presented that is intended to improve the 
sensitivity for values up to 30 GeV of the difference in mass between the two produced parti- 
cles. This broader search is motivated by the fact that an analysis performed by the ALEPH 
Collaboration, using a kinematic fit which constrains the masses of the two dijet systems to be 
equal 0, indicates that the excess events are not consistent with the hypothesis that the pro- 
duced particles have equal mass. Furthermore, compared with the ALEPH emulation analysis, 
this OPAL-specific analysis is estimated to be more sensitive to a four-jet signal of equal-mass 
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particle production at higher energies, and its efficiency is less dependent on the flavour of 
the final-state quarks. It is used in addition to the emulation of the ALEPH analysis to set 
cross-section limits as function of the dijet mass sum, and also to provide limits in the case of 
nonzero mass difference. 



2 The OPAL Detector 



A detailed description of the OPAL detector can be found elsewhere [|T0|. OPAL's nearly 
complete solid angle coverage and excellent hermeticity enable it to detect the four-jet final 
state with high efficiency. The central tracking detector consists of a two-layer silicon microstrip 
detector with polar angleQ coverage | cos^| < 0.9, immediately surrounding the beam-pipe, 
followed by a high-precision vertex drift chamber, a large-volume jet chamber and z-chambers, 
all in a uniform 0.435 T axial magnetic field. A lead-glass electromagnetic calorimeter is 
located outside the magnet coil, which, in combination with the forward calorimeter, gamma 
catcher and silicon-tungsten luminometer |]12[, complete the geometrical acceptance down to 
24 mrad from the beam direction. The silicon-tungsten luminometer serves to measure the 
integrated luminosity using small-angle Bhabha scattering events |13|. The magnet return yoke 



is instrumented with streamer tubes for hadron calorimetry and is surrounded by several layers 
of muon chambers. 



3 Data and Monte Carlo simulations 



The data used in these analyses were collected in five separate running periods. The energies 
and integrated luminosities |]T3| for the five data samples are given in Table The 130 and 
136 GeV data of 1995 and 1997 are collectively referred to in this letter as the 133 GeV data. 
The other three samples are referred to as the 161 GeV data, the 172 GeV data, and the 
183 GeV data, and are analysed separately. 

The main backgrounds for the selection of anomalous four-jet events are Z°/7* —>■ qq pro- 
duction and Standard Model four-fermion production processes. Monte Carlo samples mod- 
elling the backgrounds have been prepared using PYTHIA 5.7 fl^ for the Z^/j* qq process 
and EXCALIBUR |T^ and grc4f |T^ for the Standard Model four-fermion processes, all using 



JETSET 7.4's parton shower and hadronization models [|T^. For the generation of Standard 
Model four-fermion processes, the W mass is taken to be 80.33 GeV. Two-photon processes gen- 



erated by PYTHIA, HERWIG |8[, and PHOJET ]T9[ were used to estimate the contribution 
of these processes to the Standard Model background in the early stages of the analysis. 



The signal detection efficiencies were estimated using the HZHA generator |^ to simulate 
the production of supersymmetric Higgs bosons e"'"e~ — >■ h^A*^ — > bbbb as a model for the 



^OPAL uses a right-handed coordinate system where the +z direction is along the electron beam and where 
-t-a; points to the centre of the LEP ring. The polar angle, 9, is defined with respect to the +z direction and 
the azimuthal angle, 0, with respect to the +x direction. 
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y/ s [^ew ) 


Year 


JLat [po ) 


ioU.o 


1 one; 


Z. ( 


136.2 


1995 


2.6 


130.0 


1997 


2.6 


136.0 


1997 


3.4 


161.3 


1996 


10.0 


172.1 


1996 


10.3 


182.7 


1997 


57.1 



Table 1: Summary of the data samples, luminosity-weighted centre-of-mass energies, year of 
collection, and integrated luminosities used in these analyses. 



signal process e+e" — > XY 4 jets. Samples with decays into other quark flavours were also 
used to check for flavour dependence. All Monte Carlo samples were processed through a full 
simulation of the OPAL detector ^1 . 



4 Analysis and Results 



The main features of the signal process are four well-deflned, energetic, hadronic jets and a 
total visible event energy close to the centre-of-mass energy. The Standard Model background 
expectation changes considerably in size and composition between =133 GeV and 183 GeV. 
At 133 GeV, the main background comes from T? j^* — > qq both with or without initial-state 
radiation and accompanied by hard gluon emission. Above the threshold for e"'"e~ — > W"'"W~ 
at ^/s = 161 GeV, the background from Standard Model four-fermion processes is important 
and becomes larger with increasing y^. A procedure to reject W+W^ qqqq is implemented 
for centre-of-mass energies of 161 GeV and above. 

Events are reconstructed from charged particle tracks and energy deposits ( "clusters" ) in the 
electromagnetic and hadronic calorimeters. Tracks are required to originate from close to the 
interaction point, to have more than a minimum number of hits in the jet chamber, and to have 



a transverse momentum greater than 0.1 GeV and a total momentum less than 100 GeV p2 
Energy clusters in the electromagnetic and hadron calorimeters are required to exceed minimum 
energy thresholds. Tracks and clusters passing these quality requirements are then processed 
to reduce double- counting of energy and momentum in the event by matching charged tracks 



with calorimeter clusters. The energy-momentum flow obtained with this algorithm |2^ is used 
throughout the analysis. Energy measured in the forward detectors, covering | cos^| > 0.985, 
has not been included in the analyses presented here. 

Selection criteria emulating the ALEPH analysis ^ will flrst be described followed by the 
description of another set of requirements for event selection in an OPAL-speciflc analysis. 
Efficiencies and backgrounds for the two analyses are given followed by estimates of systematic 
errors on these quantities. 
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4.1 OPAL Emulation of the ALEPH Selection 



The procedure for selecting four-jet hadronic events and reconstructing the dijet masses is 
described below, emulating the ALEPH selection as described in Reference |Q and subsequent 
modifications and developments as described in Reference 0. The number of events retained 
after each cut in sequence is given in Table together with the expectation from Z°/7* — »• qq, 
Standard Model four-fermion and two-photon background processes, for ^/s = 133 GeV and 
for the sum over all centre-of-mass energies. Table ^ also lists the efficiencies for the reference 
h°A° Monte Carlo at ^/s = 133 GeV after each step. 



Events are required to have at least five charged tracks and seven electromagnetic calorime- 
ter clusters. The sum of the electromagnetic calorimeter cluster energies should be at least 
10% of the centre-of-mass energy, and the electromagnetic calorimeter energy is required 
to be roughly balanced along the beam direction: J2 Ei cos 9i < 0.65 Ei, where the sums 
run over measured electromagnetic calorimeter clusters. The properties of this selection 
are detailed in Reference . 



To remove events with a real Z° and large initial-state radiation (radiative return events), 
events must satisfy < K{M^is — 90 (GeV)), where is the momentum sum along 



the beam direction and Mvis = y E:^^^ — p^^^ is the total observed mass, i^' is a coefficient 
depending on ^/s. For the 133 GeV sample, K = 0.75; for the 161 GeV sample, K = 1.50; 
and for the 172 GeV and 183 GeV samples, K = 1.65. 



3. Jets are formed with the Durham jet-finding algorithm with its resolution parameter 
j/cut set to 0.008. Selected events are required to have four or more jets. For events with 
five or more jets, the jet pair with the smallest invariant mass is combined into a single 
jet and this procedure is repeated until four jets are left. 

4. The contribution of radiative return events is further reduced by requiring for each jet 
that the energy observed in the electromagnetic calorimeter, after subtracting the energy 
expected to have been deposited by the jet's charged hadrons, is less than 80% of the jet 
energy. The expected hadronic energy in the electromagnetic calorimeter is calculated 



using a track-cluster matching algorithm [23 . 



5. All jet masses are required to exceed 1.0 GeV to suppress further the contribution from 
radiative return events. 

6. The energies and momenta of each jet are rescaled imposing conservation of energy and 
momentum for the event using the beam energy constraint. The jet velocities jSi = Pi/Ei 
are held fixed in the scaling. If one or more scaling factors is negative, the event is 
rejected. The rescaled jet energies and momenta are used in the following stages of the 
selection. 

7. To suppress events involving gluon radiation, all two-jet combinations are required to 
have an invariant mass of more than 19.2% of the centre-of-mass energy. This initially 
corresponded to 25 GeV for the 133 GeV sample. 
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8. All combinations of jet pairs must have a sum of the individual jet masses {Mi + Mj) > 
10 GeV. 

9. All combinations of jet pairs must have a total charged multiplicity of at least 10. 

10. W-pair veto: For the 161 GeV data sample, a requirement is placed on the dijet mass 
sum for each of the three possible pairings of jets: M < 150 GeV for the pairing with the 
smallest AM, M < 152 GeV for the pairing with intermediate AM, and M < 156 GeV 
for the pairing with the largest AM. For the 172 and 183 GeV data, it is required that 
|M — 160 1 > 10 GeV for the pairing with the smallest mass difference if AM is less than 
15 GeV, and the same condition is applied to the pairing with the second-smallest mass 
difference if the second smallest AM is less than 30 GeV. No W-pair veto is applied to 
the 133 GeV sample. 



Cut 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Vs = 133 GeV 






















Z77* ^ qq 


3151. 


1185. 


153.2 


55.1 


50.8 


45.7 


22.5 


18.0 


14.6 


14.6 


4-Fermion 


24.5 


16.1 


4.2 


3.6 


3.1 


3.0 


1.6 


1.1 


0.7 


0.7 


77 


161.0 


8.9 


1.7 


1.7 


1.7 


1.7 


1.1 


0.0 


0.0 


0.0 


Total SM backg. 


3337. 


1211. 


159.2 


60.4 


55.6 


50.4 


25.2 


19.1 


15.3 


15.3 


Observed 


3372 


1165 


147 


51 


47 


42 


19 


16 


13 


13 


Sig. effic, ehA (%) 


99.8 


90.2 


69.0 


68.6 


67.4 


67.2 


65.4 


48.6 


45.4 


45.4 


All 






















130 184 GeV 






















Z77* ^ qq 


10456. 


4450. 


481.1 


188.7 


176.7 


155.0 


68.2 


59.7 


49.2 


40.2 


4-Fermion 


1098. 


913.0 


476.0 


428.0 


400.5 


364.8 


258.1 


236.3 


191.7 


74.2 


77 


450.7 


26.3 


2.7 


2.2 


2.1 


1.8 


1.0 


0.0 


0.0 


0.0 


Total SM backg. 


12004. 


5389. 


959.9 


618.9 


579.4 


521.6 


327.4 


295.9 


240.8 


114.3 


Observed 


12617 


5461 


984 


635 


592 


527 


328 


299 


240 


92 



Table 2: Event counts observed by OPAL at the various selection stages, with backgrounds 
estimated using PYTHIA for ZVt* qq, EXCALIBUR and grc4f for Standard Model four- 
fermion processes, and PYTHIA and PHOJET for two-photon processes. Signal efficiencies at 
133 GeV for h^A^ (see text) for M^o = Mao = 55 GeV are also listed. 



The dijet mass sum for the combination with the smallest dijet mass difference is shown in 
Figure |I|, separately for the data samples at the four different values of y/s, after the W-pair 
veto. The expected Standard Model background distribution is shown with the data for each 
case. No significant excess is observed at any of the centre-of-mass energies. 

The sensitivity of the analysis to a peak at a particular dijet mass sum depends on the 
resolution and may be affected by energy scale biases. The resolution was investigated using 
the HZHA event generator to model the process e"'"e~ — > h^A". The masses of both the h° 
and the A° are taken to be equal to 55 GeV with negligible width. The resolutions found for 
the reconstructed dijet mass sum M for the combination with the smallest AM are um = 2.0, 
2.8, 3.0, and 3.0 GeV for = 133 GeV, 161 GeV, 172 GeV, and 183 GeV, respectively. These 
resolutions do not include effects of significant non-Gaussian tails which arise from wrong jet- 
pair combinations, where the correct jet-pair combination is when each of the two jets comes 
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from the decay of the same particle. For example, at a/s = 133 GeV, 38% of the events fall into 
these tails. This definition of the resolution and tails is the same as that used in the ALEPH 
publication [T[] and the resolution values found are similar to those of ALEPH. The degradation 
of mass resolution with increasing energy arises from the scaling of the jet energies to the beam 
energy and also from the energy dependence of the detector resolution. 

Studies of the h°A° signal Monte Carlo with samples generated with input masses adding 
to 110 GeV at each centre-of-mass energy show that the reconstructed M distributions have 
peaks at masses consistent with this input value within their errors of approximately 0.5 GeV. 
Studies of events from radiative returns to the Z°, qq7, were also used to check that the Z° peak 
is well simulated in position and shape, further indicating that there is no significant bias in 
M or degradation in resolution inherent to the selection or the mass reconstruction procedure. 

For the 133 GeV signal Monte Carlo, if the reconstructed dijet mass sum for the jet-pair 
combination having the smallest AM is required to be within 2a m (4.0 GeV) of the generated 
mass sum, the efficiency obtained is 26.6%, which is 60% of the efficiency obtained without 
the requirement on M. The efficiencies and expected backgrounds both before and after the 
mass window cut are similar to those obtained by ALEPH [|l|] so that for the same integrated 
luminosity, the observed number of events can be directly compared to the number observed 
in Reference [Q. At higher centre-of-mass energies, accepting events only in a mass window of 
width ±2crM results in efficiencies of 18-21% which is 66-69% of the efficiency before the mass 
window requirement. 

To search for an excess of four-jet events with dijet mass sums near 105 GeV as motivated 
by Reference |1[], events satisfying \M — 105 GeV| < 2a m for the combination with the smallest 
AM were counted and the Standard Model backgrounds were estimated. These mass windows 
are shown in Figure |^ and the results of the searches within these mass regions are given in 
Table |^ both before and after the W-pair veto, when applicable. No significant excess is seen in 
any sample. Combining data from all centre-of-mass energies, nine events are observed while 
11.5 ± 0.4 are expected from Standard Model processes. 





Without W-Pair Veto 


With W-Pair Veto 


Data Sample 


Observed 


Expected 


Sig. eff., ehA 


Observed 


Expected 


Sig. efi:., ehA 


133 GeV 


1 


1.7 ±0.2 


26.6% 


1 


1.7 ±0.2 


26.6% 


161 GeV 


1 


1.5 ±0.1 


26.2% 





1.0 ±0.1 


18.4% 


172 GeV 


3 


2.9 ±0.1 


28.3% 





1.8 ±0.1 


20.7% 


183 GeV 


13 


11.5 ±0.4 


26.0% 


8 


7.0 ±0.3 


18.3% 


Total 


18 


17.7 ±0.5 




9 


11.5 ±0.4 





Table 3: Observed event count and expected Standard Model background for selected events 
close to 105 GeV, for the combination with the smallest AM, before and after the W-pair veto. 
No W-pair veto is applied for the 133 GeV data. The mass window is chosen to allow events that 
are within ±2aM of 105 GeV to be included, where au is the expected experimental resolution 
on M as given in the text. Signal efficiencies apply to h^A^ production with M^o = M^o = 55 
GeV. 
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Figure |^ shows the distribution of the dijet mass sum for the jet pairing with the smallest 
AM for all running periods combined, with and without the W-pair veto. The data agree 
well with the Standard Model background simulation and no excess is observed in the region 
99.0 < M < 111.0 GeV, where the width has been chosen to accommodate the resolution at 
the highest energy. A clear peak may be seen at twice the W mass in the sample before the 
W-pair veto. 

To test for a peak in the dijet mass sum distribution for arbitrary mass M and independent 
of histogram binning, the positions of the mass windows were scanned over the full range of M. 
The results are shown in Figure 0(c) for the combined data samples. The figure displays the 
event counts within windows of fixed width but whose centres are adjusted in steps of 50 MeV. 
The width of the mass window is ±4.0, ±5.6, ±6.0 GeV, and ±6.0 GeV for the 133, 161, 172, 
and 183 GeV data samples, respectively, to reflect the resolution. The contents of nearby bins 
in these scans have high statistical correlations. No significant excess is observed in the mass 
window scan at any value of the dijet mass sum. In particular, no choice of binning produces 
a peak near 105 GeV. 

To check for a possible signal in the fraction of events with wrong jet-pair combinations, 
the dijet mass sum for the jet pairing with the second-smallest AM was also considered. If 
the mass difference of a pair of objects produced together were 20 GeV, the correct jet pairing 
would yield the smallest AM for roughly half of the signal, and the second- smallest AM for 
most of the remainder. In the ALEPH analysis |jl|, including the second combination to the 
dijet mass sum distribution resulted in three additional events within the mass window with 
an additional 1.2 events expected from Standard Model processes. Figure Q shows the effect of 
adding the dijet mass sum distributions for the smallest and second- smallest AM combinations 
for different centre-of-mass energies. The distributions agree well with the Standard Model 
prediction and no peak arises when the second combination is included. 



4.2 OPAL-Specific Analysis 



In the above analysis, an emulation of the ALEPH selection criteria was applied to the OPAL 
data to test for the presence of events of the type observed by ALEPH. The selection described 
below is an OPAL-specific analysis in which the sensitivity has been maximised for detecting a 
possible signal for the process e"'"e~— >XY in the form of an excess of events with similar mass 
sums M = Mx±My in the four-jet topology. The analysis is designed to retain sensitivity even 
when the mass difference AM = |Mx — My] is as large as 30 GeV. Efficiencies and backgrounds 
are estimated for different values of M and AM. 

The cuts are designed to be as insensitive as possible to the flavours of the final state quarks. 
Although the methods employed at each of the centre-of-mass energies are similar, the optimal 
cut values in most cases depend on ^/s. 



1. The events must pass the hadronic final state requirement of cut 1 in Section 



2. The effective centre-of-mass energy after initial-state radiation, ys', calculated using the 
method described in Reference [BB], has to be at least O.STi/i. The measured visible mass. 



in 



Mvis, is required to be between i/s— 40 GeV and GeV at 133 GeV, between 100 

and 200 GeV at y/s= 161 GeV, between 110 and 210 GeV at ^/s= 172 GeV and between 
120 and 220 GeV at ^/s= 183 GeV. 

3. The charged particles and calorimeter clusters are grouped into four jets using the Durham 
algorithm [|2^. The jet resolution parameter, 1/34, at which the number of jets changes 



from three to four, is required to be larger than 0.007 at ^/s = 133 GeV, and larger 
than 0.005 at ^/s =161-183 GeV. To discriminate against poorly reconstructed events, 
a kinematic fit imposing energy and momentum conservation is required to yield a 
probability larger than 0.01. Each of the four jets is required to contain at least two tracks 
at 133 GeV and at least one track at higher energies. These kinematically constrained 
jets are used in the subsequent calculation of dijet masses. 

4. In the case of the 161-183 GeV data, the background from e'^e^ ^Z^'-f is further reduced 
by eliminating those events where one of the four jets has properties compatible with 
those of a radiative photon, namely that it has exactly one electromagnetic cluster, not 
more than two tracks (possibly from a photon conversion), and energy between 45 and 
65 GeV at y/s= 161 GeV, between 52 and 72 GeV at ■\/i= 172 GeV and between 60 and 
80 GeV at 183 GeV. 

5. The polar angle of the thrust axis, 6'thr, is required to satisfy | cos^thrl < 0.9 at 133 GeV 
and I cos^thrl < 0.8 at 161-183 GeV. 



6. To reduce background from qq events, the event shape parameter C [|23, which ranges 
between and 1 and is for a perfect 2-jet event, is required to be larger than 0.7 at 
^/s =133 GeV and larger than 0.6 at higher energies. 

7. To ensure well-separated jets for better kinematic fits, the angle between any two jets is 
required to exceed 0.8 radians for 161-183 GeV data. 

8. Above the W"^W~ threshold, explicit vetoes against the process e"'"e~^W''"W~ are ap- 
plied. 

At ^/s =161 GeV, the two bosons are produced with only a small boost. The two 
jets having the largest opening angle are assigned to one of the bosons and the two 
remaining jets to the other. An event is rejected if both jet pairs have an invariant mass 
between 75 GeV and 90 GeV. 

At y/s =172 and 183 GeV, a more sophisticated veto is applied. The four jets are 
combined into pairs, and for all three combinations the event is refitted constraining the 
total energy to ^/s and the total momentum to zero, and also constraining the masses 
of the two jet pairs to be equal (five constraints). From the three combinations, the 
one yielding the largest fit probability is considered. If the jet pair mass from the fit 
exceeds 75 GeV and the fit probability is at least 0.01, the event is rejected. 

9. To achieve good sensitivity for all AM less than 30 GeV, we use two separate mass 
selections, one relevant for unequal masses and one for equal masses. In both selections, 
when searching for a signal with a hypothetical sum of masses, Mq, the range Mq ± 2a m 
is used, where ctm is 2.0 GeV at ^/s =133 GeV and 3.0 GeV for higher energies. For 
unequal masses (AM > 5 GeV), the event is selected if either the jet association with 
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the smallest mass difference or the one with the second smallest mass difference has a 
mass sum M in the range Mq ± 2a m- For nearly equal masses (AM < 5 GeV), better 
sensitivity is obtained when considering only the jet association with the smallest mass 
difference. The resolution ctm varies only slowly with M and AM. 

Table ^ presents the number of observed events and the Standard Model expectations before 
and after the W"'"W~ veto (cut 8). The numbers of observed events are consistent with the 
Standard Model expectations at all centre-of-mass energies both before and after the mass 
selection. 





Without W-Pair Veto 


With W-Pair Veto 


After Mass Selection 


Data Sample 


Observed 


Expected 


Observed 


Expected 


Observed 


Expected 


133 GeV 


18 


17.0 ±0.6 


18 


17.0 ±0.6 


4 


3.1 ±0.3 


161 GeV 


11 


15.8 ±0.3 


8 


13.6 ±0.3 


2 


2.7±0.1 


172 GeV 


36 


33.8 ±0.3 


21 


16.2 ±0.2 


4 


2.9 ±0.1 


183 GeV 


190 


210.1 ± 1.2 


70 


81.6 ±0.8 


6 


8.9 ±0.3 


Total 


255 


276.7 ± 1.4 


117 


128.4 ± 1.1 


16 


17.6 ±0.4 



Table 4: Number of observed and expected Standard Model background events before and 
after the W-pair veto and after adding the mass selection (cut 9) centred at 105 GeV for the 
smallest AM combination. The quoted errors are statistical. No W-Pair veto has been applied 
to the 133 GeV data. 



Table ^ shows the signal efficiencies for various combinations of (Mx,My) together with the 
predicted background and the numbers of observed events after all cuts. 

Figure § shows the distributions of M for the jet associations with the smallest AM, and 
for the jet association with the smallest and second-smallest AM, summed over all centre-of- 
mass energies. Globally, the distributions show consistency between the data and the Standard 
Model background prediction. In particular, there is no excess in the vicinity of M ^ 105 GeV. 
The overlap of the OPAL-specific analysis and the OPAL emulation of the ALEPH selection 
has been evaluated in a typical Monte Carlo four-jet signal sample at ^/s = 133 GeV with 
Mho = Mao = 55 GeV. In this sample, 59% of the events selected by the OPAL emulation of 
the ALEPH analysis are also selected by the OPAL-specific analysis. 

4.3 Systematic Errors 

At ^/s = 133 GeV, since the efficiencies and expected backgrounds for the OPAL emulation of 
the ALEPH signal are similar to those obtained by ALEPH it is not necessary to consider 
systematic effects in detail if only numbers of observed events are compared. However, to cal- 
culate limits on cross-sections, systematic errors on efficiencies and backgrounds are estimated. 

To emulate the ALEPH analysis, the charged multiplicity requirement (cut 9, Section 4.1) 
on all combinations of two jets was necessary. Since the aim of the emulation analysis is to 
compare directly the OPAL result with the ALEPH observation, and ALEPH's cross-section 
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estimate was made assuming a model of four b-jets, we also assume this model. Varying the 
mean charged multiplicity in b-hadron decays by its measurement uncertainty results in an 
estimated systematic error of 12% on signal detection efficiencies due to this effect. Including 
additional uncertainties in the modelling of the cut variables, energy scales, mass resolutions 
and limited Monte Carlo statistics results in an estimated total systematic error of 13%. 





133 GeV 


161 GeV 


(Mx,My) 


Effic. 


Backgd. 


Data 


Effic. 


Backgd. 


Data 


(GeV) 


(%) 






(%) 






(50,50) 


30.6 ± 1.5 


3.2 ±0.3 


3 


38.0 ±2.2 


2.4 ±0.1 


3 


(40,60) 


29.6 ±2.0 


5.6 ±0.4 


6 


35.8 ±2.1 


4.2 ±0.2 


3 


(55,55) 


29.0 ±2.0 


3.1 ±0.3 


4 


37.4 ± 1.0 


3.0 ±0.1 


1 


(50,60) 


38.4 ± 1.5 


5.6 ±0.4 


7 


42.0 ±2.2 


5.0 ±0.2 


1 


(40,70) 


23.0 ± 1.9 


5.6 ±0.4 


7 


34.8 ±2.1 


5.0 ±0.2 


1 


(60,60) 


26.2 ± 1.4 


3.2 ±0.3 


3 


39.0 ±2.2 


2.8 ±0.1 





(50,70) 


30.2 ±2.1 


5.6 ±0.4 


4 


41.2 ±2.2 


4.8 ±0.2 


1 


(60,70) 


24.2 ± 1.9 


2.2 ±0.2 


4 


34.2 ±2.1 


4.7 ±0.2 


2 


(50,80) 


13.1 ± 1.1 


2.2 ±0.2 


4 


32.8 ±2.1 


4.7 ±0.2 


2 


(70,70) 








26.6 ±2.0 


2.0 ±0.1 


2 


(60,80) 








33.8 ±2.1 


4.4 ±0.2 


3 





172 GeV 


183 GeV 


(Mx,My) 


Effic. 


Backgd. 


Data 


Effic. 


Backgd. 


Data 


(GeV) 


(%) 






(%) 






(50,50) 


31.0 ± 1.5 


2.5 ±0.1 


4 


19.4 ± 1.8 


4.8 ±0.2 


4 


(40,60) 


24.6 ± 1.9 


3.2 ±0.1 


5 


16.3 ± 1.6 


5.6 ±0.2 


6 


(55,55) 


34.4 ± 1.0 


3.6 ±0.2 


5 


31.6 ±2.1 


12.4 ±0.3 


9 


(50,60) 


32.6 ± 1.5 


5.1 ±0.2 


6 


22.2 ± 1.9 


16.2 ±0.3 


12 


(40,70) 


23.0 ± 1.3 


5.1 ±0.2 


6 


18.0 ± 1.7 


16.2 ±0.3 


12 


(60,60) 


33.0 ± 1.5 


3.7 ±0.2 


9 


32.6 ±2.1 


17.1 ±0.4 


20 


(50,70) 


34.4 ±2.1 


5.5 ±0.2 


11 


23.9 ± 1.9 


23.5 ±0.4 


30 


(60,70) 


33.6 ± 1.5 


5.5 ±0.2 


6 


32.2 ±2.1 


28.2 ±0.4 


24 


(50,80) 


28.8 ± 1.4 


5.5 ±0.2 


6 


23.2 ± 1.9 


28.2 ±0.4 


24 


(70,70) 


29.9 ± 1.4 


3.3 ±0.1 


1 


31.2 ±2.1 


20.7 ±0.4 


20 


(60,80) 


31.0 ±2.1 


6.0 ±0.2 


4 


27.7 ±2.0 


32.3 ±0.5 


28 



Table 5: Signal detection efEciencies, numbers of expected background events and number of 
observed data events, for various mass combinations in the OPAL-specific analysis, after the 
mass selection, cut 9 of section 4.2. The quoted errors are statistical. 

In the OPAL-specific analysis, the signal detection efficiencies are subject to a systematic 
error of 9%, which includes an allowance for the final state to contain any composition of 
quark flavours and uncertainties in modelling heavy hadron decays, the uncertainty on the 
simulation of the decay with regards to fragmentation and hadronization, the modelling of the 
cut variables, and the limited Monte Carlo statistics. 
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The total relative uncertainty on the residual background is 20% for the OPAL emulation 
of the ALEPH analysis, and 13% for the OPAL-specific analysis. These errors include the 
uncertainty on the modelling of the hadronization process, on the prediction of the four-jet 
rate, W-pair cross-section, and the modelling of the cut variables. The error due to the limited 
Monte Carlo statistics is added in quadrature to this uncertainty. The systematic errors on the 
luminosity measurements range from 0.5% to 1.6%. 



5 Cross-Section Upper Limits 

In the OPAL emulation of the ALEPH analysis, the number of observed events can be compared 
directly to the ALEPH observation ||l| because both the observed background rate and the 
estimated efficiency are nearly identical to those obtained by ALEPH. From the number of 
observed and expected events in the dijet mass sum window of 105±4 GeV at y/s ^ 133 GeV, we 
set a 95% confidence level (CL) upper limit of 2.1 events that could be attributed to additional 
cross-section from new physics when scaled to the integrated luminosity of the 1995 ALEPH 
result. This can be compared to ALEPH's observation in 1995 of nine events with a Standard 
Model expectation of 0.8 events. To calculate the probability that the OPAL observation is 
consistent with the ALEPH observation in the presence of a possible signal, the product is 
formed of the Poisson probability pi that at least nine events were observed in ALEPH and p2 
that no more than one event was observed in OPAL, given the Standard Model backgrounds 
and assuming the presence of a signal scaled by the integrated luminosity. The probability of an 
outcome no more likely than that observed in the data, i.e., the sum of Poisson probabilities of 
possible outcomes less than or equal to piP2, is found to be 2.6x10"^, where the hypothesized 
signal cross-section has been chosen to maximize this probability. 

Assuming production of new particles X and Y subsequently decaying to a final state of four 
b-jets to determine efficiencies, cross-section upper limits are set using the OPAL emulation 
of the ALEPH analysis. From the number of observed and expected events in the dijet mass 
sum window as above at ~ 133 GeV, a 95% confidence level (CL) upper limit of 1.4 pb is 
determined for the production cross-section at the dijet mass sum of 105 GeV. Systematic un- 
certainties of efficiency, background and luminosity are taken into account using the procedure 
outlined in Reference |l29| . 

To combine data from different centre-of-mass energies, we consider two different functions 
for the energy dependence of the cross-section of a hypothetical signal. It is first assumed 
that the cross-section varies as P{3 — l3^)/s, typical for pair-production of spin-1/2 particles, 
where f3 is taken as the average velocity of the particles in the laboratory frame Taking 
from Table ^ the total number of observed and expected events in resolution-dependent mass 
windows around 105 GeV, an upper limit on the production cross-section at 133 GeV of 0.58 pb 
at 95% CL is found. Secondly, under the hypothesis that the signal cross-section varies as 
typical for the production of scalar particles , the upper limit on the cross-section at 133 GeV 



is computed to be 0.31 pb at 95% CL. These limits can be compared to ALEPH's estimated 
cross-section of 3.1 ± 1.7 pb pi] from their total number of excess events. 



U 



The OPAL-specific analysis described in section 4.2 is used to obtain upper limits for the 
cross-section of a possible signal process e"'"e~— s>XY— > four jets, in the presence of background 
from Standard Model processes, using Poisson statistics and incorporating systematic uncer- 
tainties as described in Reference p9|. The process e"'"e~— i>h°A° was used to model the sig- 



nal detection efficiencies. The resulting 95% CL upper limits, as function of the mass sum 
M(= Mx -|- My), are shown in Figure ^, for AM close to zero and AM = 30 GeV. A mass 
window of M ± 2(Tm is scanned across the distribution of the dijet mass sum in small steps. 
To account for a possible discrepancy between the mass scale of the data and the Monte Carlo 
in a conservative manner, the mass window is displaced by ±0.5 GeV at each scan point. The 
largest data count in any of the three windows including the nominal one and the smallest back- 
ground estimation in any of the three windows are used to compute the limit. When results 
at different centre-of-mass energies are combined, the hypothetical production cross-section is 
assumed to vary as The cross-section limits are presented separately for the 133 GeV 

data, and for all data (130-184 GeV) combined. Limits on the cross-section from the combined 
data sample are computed both at y/s = 133 GeV and at y/s = 183 GeV. These results are 
independent of the flavour of the quarks from the decay of the hypothesized particles and are 
valid for X and Y being scalars produced predominantly by an s-channel process. 



6 Conclusions 



Following the ALEPH observation of a large excess of four-jet events with dijet mass sums 
around 105 GeV at ^/s ~ 133 GeV IQ, a careful emulation of the ALEPH analysis has been 
performed using OPAL data collected from e+e~ collisions at centre-of-mass energies between 
130 and 184 GeV. The process e+e~ —>■ h°A° was used to estimate the signal detection efficien- 
cies. The estimated sensitivity, mass resolution, efficiency, and estimated backgrounds in this 
analysis were similar to that of the ALEPH analysis. No significant excess of four-jet events 
with dijet mass sums in the region close to 105 GeV, or any other region between 60 and 160 
GeV, has been observed in any of the data samples separately or combined, and our obser- 
vations are consistent with Standard Model predictions. The same conclusions are reached 
when an OPAL-specific analysis is employed. Limits for the cross-section of a hypothetical 
process e"'"e~^XY— >■ four jets are given as a function of the dijet mass sum M and the mass 
difference AM. The 95% confidence upper limits obtained in both analyses for dijet mass 
sums near 105 GeV are below the excess reported in 1995 by ALEPH |jl| to a high degree of 
confidence. ALEPH has recently analysed new data at centre-of-mass energies between 
130 and 184 GeV and do not confirm the previously reported excess. 
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Figure 1 : The dijet mass sum for the combination with the smallest AM after the W-pair veto 
in opal's emulation of the ALEPH analysis shown separately for the different centre-of-mass 
energies. Data are shown by the points and Standard Model backgrounds by the histograms. 
The hatched component of the background histograms denotes Standard Model four-fcrmion 
processes, while the unhatched component denotes T? j^* qq. The mass windows containing 
the region of interest are indicated by the arrows. The clashed histogram in (a ) illustrates a 
signal (plus background) that could be expected due to h°A° with both decaying to pairs of 
b-quark jets, and M^fi — M^o — 52.5 GeV, nem^alized to the excess observed by ALEPH at 
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Figure 2: The dijet mass sum in OPAL's emulation of the ALEPH analysis for the combined 

130-184 GcV samples. Plots (a) and (h) show the distribution of the dijet mass sum before 
and after the W-pair veto, respectively. Data are shown by the points and Standard Model 
backgrounds by the histograms. The hatched component of the background histograms denotes 
Standard Model four-fermion processes and the unhatched component denotes Z^/Y — QQ- 
The mass window around 105 GeV whose width accommodates the resolution at = 183 GeV, 
is shown with the arrows. Plot (c) shows the sliding mass window scan for the same analysis 
after the W-pair veto. The hatched histograms show the total number of data events, and 
the solid line shows the Standard Model expectation; the line width indicates the Monte Carlo 
statistical error. An arrow is drawn at 105 GeV. 
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Figure 3: The dijet mass sum in OPAL's emulation of the ALEPH analysis for both the 
combination with the smallest AM and the combination with the second-smallest AM for 
(a) the 133 GeV data sample and (b) the combined 130-184 GeV samples. The points and 
histograms are as in Figure [J. 
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Figure 4: Distributions of M in the OPAL-speciGc analysis for the combined 130-184 GeV sam- 
ples after all selection requirements except the mass selection (cut 9), (a) for the jet combination 
with the smallest AM and (b ) for the jet combinations with the smallest and second-smallest 
AM. The points and histograms are as in Figure 
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Figure 5: The 95% CL upper limits obtained with the OPAL-specihc analysis on the production 

cross-section of a possible signal as a function of M for AM close to (solid lines) and for AM 
< 30 GeV (dashed lines). Plot (a) shows the limits computed using the data collected at 
^/s^l33 GeV; plot (b) shows the limits using the combined data from ^/s = 130-184 GeV 
assuming a cross-section that varies as P'^/s scaled to y/s = 133 GeV (lines that end near 
M=130 GeV) and ^/s = 183 GeV (lines that extend to larger M). 
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